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Abstract 

The response of the hquid scintillator array Nordball to neutrons in the energy 
range 1.5 < T„ < 10 MeV has been measured by time of flight using a ^^^Cf fission 
source. Fission fragments were detected by means of a thin-film plastic scintillator. 
The measured differential and integral neutron detection efficiencies agree well with 
predictions of a Monte Carlo simulation of the detector which models geometry ac- 
curately and incorporates the measured, non-linear proton light output as a function 
of energy. The ability of the model to provide systematic corrections to photoneutron 
cross sections, measured by Nordball at low energy, is tested in a measurement of the 
two-body deuteron photodisintegration cross section in the range = 14 — 18 MeV. 
After correction the present (7, n) p measurements agree well with a published 
evaluation of the large body of (7, p) n data. 
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1 Introduction 



Recently photoneutron production cross sections on a range of nuclei [1] have 
been measured in the energy region 11-30 MeV at the high-duty- factor, tagged- 
photon facility of MAX-lab [2] in Lund, Sweden. The purpose is to test cal- 
culations of neutron dose received during the course of bremsstrahlung ra- 
diotherapy [3]. In these time-of- flight (TOF) experiments, at relatively low 
kinetic energy (T„), the neutron signal is obscured by an accelerator-induced 
room background and a high rate of random-coincidence events produced by 
untagged photons incident on the experimental target. The bulk of this back- 
ground arises from detected photons, so that good particle identification is 
required to access T„ ~ 1 MeV, where both the photoproduction cross section 
and the biological effectiveness of neutrons are at a maximum. The MAX-lab 
liquid scintillator array Nordball [4], with good n/7 pulse shape discrimina- 
tion (PSD) properties and ~ 1 ns (FWHM) time resolution, was used for 
these measurements. This paper presents Nordball calibration procedures and 
a computer model used to evaluate systematic effects which distort the mea- 
sured neutron yield. 

At low energies neutron attenuation and multiple scattering effects on mea- 
sured neutron yields are large. Thus a simulation of the experiment based 
on GEANT-3 [5], which models neutron interactions in all materials in the 
vicinity of the detector array, has been developed to correct for these effects. 
The non-linear pulse-height response of the scintillators to low energy recoil 
protons has a critical bearing on these calculations, especially close to thresh- 
old, and so this was measured with a white ^^^Cf neutron source. This is a 
calibration standard since the neutron yield from spontaneous fission (decay 
branching ratio 3.1%) is high and both the absolute numbers and energy spec- 
trum of prompt neutrons per fission are well known [6,7]. Neutron detector 
response may be obtained in a TOF measurement if fission events are tagged 
by detection of at least one heavy fragment, which provides a time reference. 
Although fission chambers [8] are possibly the optimal detectors for fragment 
detection, a simple thin-film organic scintillator presented a viable, available 
alternative at the time of the measurement. This type of detector has several 
desirable properties: 

(1) Insensitivity to 7-rays or neutrons as the plastic is extremely thin. 

(2) Fast response with a similar time resolution to Nordball, which also gives 
high counting rate capability. 

(3) Fast, cheap production. 

Here we compare the measured and simulated fission-neutron response of the 
Nordball array and describe a test of the simulation by measurement of the 
well known (7, n) p cross section. Here the simulations have been used to 
compute the large attenuation and multiple scattering corrections, as well as 
the neutron detection efficiency. Sec. 2 gives an overview of the experimental 
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setup, the measured Nordball response is presented in Sec. 3, the Monte Carlo 
(MC) simulation is described in Sec. 4, measurements and simulations are 
compared in Sec. 5, and a short summary is given in Sec. 6. 

2 The time-of-flight experiment 

2.1 The Nordball array 

The Nordball detector (Fig. 1) consists of 16 liquid scintillators, type Bicron 
BC-501, with PSD capability. Ten detectors are of hexagonal cross section and 
six detectors of pentagonal cross section. Their respective volumes are 3.3 1 and 
2.6 1 at a common thickness of 16 cm. The liquid is contained within a 2 mm 
thick, stainless steel cannister, coated on its inner surface with Ti02 reflective 
paint. This is connected to a 5 in. XP2041 photomultiplier tube (PMT) via 
a Pyrex glass window. A cylinder of /i-metal shields the PMT from stray 
magnetic fields and an outer plastic housing encases tube and voltage-divider 
circuit, which provides a negative anode signal. 

All detectors were mounted on aluminium frames and placed on a 32 cm thick 
layer of borated paraffin, supported by an iron table. The configuration of 
Fig. 1 consisted of five detector columns, positioned in 15° steps, at a distance 
of 150 cm from the central position where the experimental target or the 
fission detector was located. The bigger hexagonal detectors were placed in 
the two bottom rows, five pentagonal detectors in the third, and the last one 
on top of the central column. Paraffin towers to both sides of the iron table 
partially shielded the array against regions of strong neutron background. 
For the deuteron photodisintegration experiment a 40 x 40 x 1.8 cm plastic 
scintillator sheet was inserted between the target and Nordball to identify 
charged particle events. 

2.2 Fission detector 

The thin-film scintillator was prepared following Ref. [9]. A solution of plastic 
scintillator in xylene was spread uniformly over the glass window of a 2 in. 
PMT (Philips XP2262B) which was placed upright in a vacuum chamber. 
Slow evacuation of the chamber causes solvents to evaporate and a thin layer 
of scintillating material remains on the glass. Its thickness was tuned, by vary- 
ing the amount of plastic in the solution, to optimise discrimination between 
fission fragments and the 30 times more numerous a particles. 

The fission source [10] consists of a platinum-clad, nickel disk, on to which 
Cf203 was electro-deposited. The active area (0.2 cm^) is covered by a 50 /xg/cm^ 
layer of gold allowing the passage of fission fragments with relatively small en- 
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Figure 1. Layout of Nordball array for ^^^Cf fission-source and (7, n) p measure- 
ments. 

ergy loss. The source was placed on the centre of the scintillating film and 
the PMT was sealed with a plastic cap, lined with aluminised mylar foil, to 
exclude external light and provide some reflection of scintillation light. The 
flssion detector was positioned so that the source sat at the target-centre po- 
sition (Fig. 1) at the same height as the middle of the second Nordball row 
and the photon beam axis. 

The axis of the PMT was offset ~ 10° from vertical. Apart from the bottom 
row of the array, this offset avoided fission neutrons having to pass through 
the glass of the PMT on the way to the Nordball detectors. Placing the PMT 
horizontally would have displaced the loose source from the active area of the 
PMT face. As a consequence of attenuation in glass, the detected neutron yield 
in the bottom detectors was found to be about 15% smaller at low energies 
where the interaction cross section is highest. These losses are consistent with 
the predictions of a computer model of the experiment (Sec. 4.2) which ap- 
proximated the PMT by a 2 mm thick glass cylinder and calculated for tilting 
angles in the range of 5-15 degrees. However the angle was not sufficiently 
well determined to enable quantitative comparisons and hence the bottom de- 
tectors were excluded from the analysis of differential and absolute neutron 
detection efficiencies presented in Sec. 5. 

No suppression of bottom-row neutron yield was observed in the deuteron 
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photodisintegration measurement, which used a cylindrical D2O target, sup- 
porting the assumption of increased neutron attenuation in glass during the 
^^^Cf measurement. 

Fig. 2 shows the fragment pulse-height distribution obtained from the ^^^Cf 
source, which produced a counting rate of ~ 3 kHz. With ~ 27i acceptance for 
charged particle detection, the effective thickness of the scintillator depended 
on emission angle and not all 6.11 MeV a particles could be separated from 
heavy fragments by pulse height. A small fraction of a particles had to be 
accepted in the trigger, resulting in the sharp spike observed around channel 
250. These events are completely uncorrelated with any detected neutron, 
confirming that they are not produced by fission fragments. The two broad 
distributions represent interactions of heavy and light fission fragments, with ~ 
10% of the former lost, judging from the width of the peak. An estimate of the 
resulting distortion of the neutron spectrum was based on Ref. [6] where the 
shapes of the neutron energy distributions and the mean numbers of emitted 
neutrons per fission event for light (1.97) and heavy (1.70) fragments are 
presented. The distortion in the neutron spectrum is below 2%, and the mean 
number of neutrons V is reduced by 5.4%. No significant deviations from this 
simple estimate were observed when measured neutron spectra, corresponding 
to different regions in the fission pulse-height distribution, were compared. 
However rather than attempt a more sophisticated correction we have factored 
these distortion effects into the systematic uncertainties associated with the 
evaluation of the ^^^Cf fission neutron yield [7]. 

2. 3 Electronics 

The layout of the electronics is shown in Fig. 3 and the data acquisition and 
analysis system ACQU is described in Ref. [11]. Analogue to digital converters 
(ADC) and scalers were read out by a VMEbus single board computer via the 
CAMAC bus and a VMEbus CAMAC parallel branch driver (CBD). During 
readout the event latch prevented further triggers until reset at the end of the 
event. 

The experimental trigger for the ^^^Cf measurement was made by the constant- 
fraction discriminator (CFD) of the fission detector. The number of triggers 
was recorded by a scaler for normalisation purposes and the fission pulse height 
was recorded in a charge-integrating ADC (QDC). Gate and start signals for 
the various ADCs were derived from the trigger. The anode signal of each 
Nordball scintillator was split 3 ways and fed to: 

(1) a QDC to record the pulse amplitude. 

(2) a CFD connected to a time to digital converter (TDC) to record the time 
of fiight to the scintillator. 

(3) specialist hardware developed for fast pulse-shape analysis [12], which 
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Pulse height (Channel) 

Figure 2. ^^^Cf pulse-height spectrum obtained with the thin-film fission detector. 
The narrow peak results from those alpha particles remaining above pulse-height 
threshold. The two broad distributions correspond to heavy and light fission frag- 
ments. 



produces a "pulse-shape" output, recorded in a voltage to digital converter 
(VDC). 

The electronics for the (7,n) p measurement were very similar, except that 
the trigger was made by the logical OR of the Nordball CFD outputs. In 
this case the logic outputs from the pulse-shape analysers, flagging a neutron, 
were used to gate the interrupt signal to the CPU. In addition the signals from 
each of the 64 focal-plane detectors of the photon-tagging spectrometer were 
recorded in TDCs and scalers. 
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Figure 3. Block diagram of the experimental electronics (Sec. 2.3). 
Nordball Performance 



3.1 Pulse-height calibration 



The light output from the Nordball detectors was calibrated with the 7-ray 
sources ®°Co, ^^^Th, and Pu/Be. Following Ref [13], the edges of the resulting 
Compton electron spectra were initially assumed to be the channels containing 
90% of the counts in the peaks in the Compton end-point pulse-height distri- 
butions. However this procedure may be over simplistic as the Compton edge 
position depends on the resolution of the detection system under consideration. 
In Ref [14] an empirical formula was determined to estimate the resolution 
from measured pulse-height spectra and tables of resolution-dependent Comp- 
ton edge positions, relative to the channel of maximum counts, were produced. 
These tables were checked using the present GEANT-3 simulation (Sec. 4), 
which reproduced the result for both Nordball and the detector geometries of 
Ref [14]. However the simple approximation of Ref [13] does not deviate from 
Ref [14] by more than 3%. 

The pulse-height resolution (FWHM) of the Nordball array was parameterised 
[14], in terms of scintillation light output L in units of electron energy equiv- 
alent (MeVee), as: 

^ = ^Ja^ + pyi + ^yi^ (1) 



Small non-linearities in the low energy electron response [15] were neglected. 
Since the individual detector responses were very similar, they were summed 
before least-squares fitting to obtain the average response of the Nordball 
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array. Parameter values of a = 3.2 ±10 % and P = 10.9 ±5 % were obtained 
with 7 fixed at zero and this resolution function was implemented in the MC 
simulation. The obtained values are similar to those of Ref. [16,17,18]. 

3.2 Particle identification 

Pulse shape discrimination (Fig. 4) was applied to separate detected neutrons 
and photons. To quantify neutron and photon regions, two-Gaussian fits were 
made to pulse shape projections generated for 50 keVee-wide pulse-height bins, 
yielding peak positions and widths for each distribution. Intervals of ±3cr 
around the peak positions defined the neutron or photon region limits at a 
given pulse height and straight-line fits to these limits defined the boundaries 
of the neutron/photon regions of the pulse shape plot. Neutron events below 
the dashed line of Fig. 4 have some photon contamination at pulse heights 
below ~ 1 MeVee. 

3.3 The neutron pulse-height response function 

The non-linear scintillation response to recoil protons is vital input to a neu- 
tron efficiency simulation and must be measured to avoid large systematic 
errors close to detection threshold. Using TOF information, neutron events 
were divided into 625 keV-wide kinetic-energy bins in the range 1.25 < T„ < 
10 MeV. The resulting pulse-height spectra for each detector were summed to 
produce the overall response of the array. For each value of T„ the maximum 
proton recoil energy was determined from the half-height position Ci/2{Tn), 
corrected by a multiplicative factor C(T„) to account for multiple neutron 
interactions in the scintillators. 

Numerical values for C(T„) were derived from the MC simulation (Sec. 4.1), 
which mimics the 625 keV neutron energy bins used in the data analysis. 
A pulse-height distribution for 5 MeV neutrons is shown in Fig. 5. The half- 
height position Ci/2 at 1.85 MeVee underestimates the expected light output of 
2.03 MeVee for maximum-energy recoil protons, calculated from the scintillator 
response function. This leads to a correction factor C(T„) = 1.10 at 5 MeV, 
while for 2 MeV neutrons C(T„) = 1.22. A 2nd-order polynomial fit to the 
MC calculations in the range T„ = 1.25 — 10 MeV yielded: 

C{Tn) = 1.335 - 0.067 ■ T„ + 0.004 ■ (2) 

The proton response functions obtained with (filled circles) and without (empty 
circles) the correction factors are shown in Fig. 6. The symbol size is equivalent 
to the 50 keVee uncertainty of the half-height position C1/2, estimated from 
variations in the spline functions which fitted the edge regions of the pulse- 
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Pulse height (MeVee) 

Figure 4. Neutron-photon discrimination plot for a single Nordball detector, showing 
the pulse-shape information as a function of pulse height. The dashed/full lines 
represent the boundaries of the neutron/photon regions. 

height spectra. The response function of Ref. [19] (dashed line) describes the 
energy dependence of the data quite well, after multiplication by a factor of 
0.9 (full line) which was determined from a least-squares fit. Note that the first 
data point at 1.25 MeV is too high because the hardware threshold of 0.25 
MeVee cuts into the recoil proton edge distribution (Sec. 5.1). The response 
of individual detectors showed no significant deviations and a common light 
output function for the Nordball array was implemented in the MC simulation. 

4 The GEANT-3 based Monte Carlo simulation 

Measurements of low energy neutron photoproduction require correction for 
sizable systematic effects on the neutron acceptance, which has led to the 
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Figure 5. The simulated pulse-height distribution for 5 MeV neutrons. 

development of a GEANT-3 [5] based MC simulation, into which neutron 
tracking methods from the STANTON code [20,19] have been incorporated. 
This coupled approach allows the use of the GEANT-3 neutron cross section 
data base (package GCALOR [21] based on ENDF/B-6 [22]) at T„ <20 MeV 
and STANTON procedures for 20 < T„ < 500 MeV, where the latter provides 
a more detailed and accurate representation of neutron interactions compared 
to standard GEANT-3. 

The flexible geometry framework of the MC simulation allowed comparison of 
predictions with a broad range of published neutron response measurements. 



4.1 The modified STANTON code 



A modified version of STANTON, augmented in the following ways, was used 
to model the Nordball response in the ^^^Cf measurement: 
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Proton energy (MeV) 

Figure 6. The pulse-height response function of the Nordball array: full/empty circles 
represent the half-height positions with/without correction for multiple interactions. 
The dashed line represents the response function of Ref. [19], the full black line the 
same function multiplied by a factor of 0.9. 

(1) The hexagonal Nordball geometry was added. Pentagonal elements were 
treated as hexagonal, but scaled down in cross-sectional area to give the 
correct pentagonal-element volume 

(2) Isotropic neutron emission was added to the existing event generator for 
the ^^^Cf source. 

(3) The pulse-height resolution function of Sec. 3.1 was incorporated. 

(4) The scintillator response to proton interactions, as described in Sec. 3.3, 
was incorporated. 

(5) Results were stored as CERN HBOOK ntuples for later analysis. 

The core of the simulation, including the neutron tracking routines of Ref. 
[19], remained unchanged. After benchmark comparisons which reproduced a 
variety of published neutron efficiency measurements [23,24,25], the MC code 
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was applied to Nordball (Sec. 5.1,5.2). 



4.2 The coupled GEANT3-STANT0N code 

The geometry and materials of a particular experiment (Fig. 1) can be de- 
scribed in fine detail using GEANT-3 but, designed primarily for high energy 
physics, this fails to model low energy neutron interactions accurately. How- 
ever, with the addition of the GCALOR package the GEANT-3 description 
of Tn < 20 MeV neutron interactions is quite accurate due to the use of re- 
cent Evaluated Nuclear Data Files. These files contain neutron partial cross 
sections, angular distributions and secondary energy distributions evaluated 
on a grid of discrete energy and angle points, for elements over the full range 
of the periodic table. This is necessary as few-MeV neutron cross sections are 
large and strongly dependent on energy and target mass. Any realistic model 
requires consideration not only of neutron interactions within the detector 
volume, but also of attenuation and multiple scattering effects in all bulk ma- 
terials in the locality of the detector. At higher energy such effects become 
less important, but none the less it remains important to model interactions 
in the detector volume as accurately as possible. 

The present implementation of GEANT-3 incorporates STANTON methods 
as a subroutine which tracks neutrons inside the organic-scintillator detec- 
tor volumes, while GCALOR is used for interactions outside the detector 
volume. The geometry of all components considered in the simulation, and 
hence calculation of volume boundaries, is controlled entirely by the GEANT- 
3 tracking system. Original STANTON geometry functions were replaced by 
the GEANT-3 subroutine gmedia, which may be tailored to handle pentagonal 
and hexagonal shaped detectors in a straightforward way. 

A test of the simulation is described in Sec. 5.3, where it has been used to 
calculate systematic corrections to a measurement of the well known, two-body 
deuteron photodisintegration cross section. 



5 The detection efficiency from the ^^^Cf measurement 

5.1 Differential detection efficiency 

Normalised pulse-height spectra were sorted into incident neutron energy bins 
of 0.5 MeV width and compared to the MC predictions. A PSD cut was ap- 
plied to select neutrons and events with unreasonably large recoil proton pulse 
height, compared to that expected from the incident neutron energy and mea- 
sured proton response, were rejected. 
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The resulting TOF spectrum for a hardware threshold of ~ 0.25 MeVee, shown 
in Fig. 7, is virtually free of random coincidences which would fill in the 20- 
30 ns TOF region between the photon peak (some photons beat the PSD cut) 
and the neutron distribution. A low energy tail due to neutrons that scattered 
en route to the detectors (and hence had a longer flight path) is present at flight 
times of 100 ns and longer. These events were suppressed when the detection 
threshold was raised, and above 0.7 MeVee they were virtually eliminated as 
discussed in Sec. 3.2. 




Time-of-flight (ns) 

Figure 7. Time of flight spectrum for a single Nordball detector. The width of the 
photon peak is 1.7 ns (FWHM). 

Assuming isotropic fission neutron emission and neglecting absorption effects 
in air and detector canisters, which were estimated to be below 3%, the dif- 
ferential efficiency is given by: 

dL^ V- Nf,,, ■ ■ AT„ • AQ/Att ^ ^ 



13 



where ^jf-{Tn) is the pulse-height spectrum for a neutron energy bin of width 
AT„ centered at T„, V = 3.76 (5.4%) is the average number of prompt neutrons 
per fission event, Nfi^g is the number of fission fragments counted, fn{Tn) 
(3.5%) is the energy distribution of prompt neutrons and AQ (0.4%) is the 
solid angle subtended by a detector. Relative uncertainties are given in the 
brackets. The main sources of uncertainty result from distortions of V and 
fn{Tn) due to not recording the full fission- fragment, pulse-height spectrum 
(Sec. 2.2) and the fniTn) uncertainty contains a further factor due to the 
uncertainty in T„. 

The agreement between data and simulation is generally good, as shown in 
Fig. 8. Since the statistical uncertainty increases rapidly as neutron energy 
increases, the spectra were binned a factor of 2 more coarsely above 5 MeV. 
The hardware threshold cuts noticeably into the pulse-height distribution for 
values below 0.3 MeVee. This (voltage) threshold is not sharp due to pulse 
shape dependence and hence difficult to model accurately. Thus 0.3 MeVee 
represented a minimum software threshold which could be used for quanti- 
tative analysis. At low incident energies multiple scattering leads to a more 
peak-like shape of the spectrum which flattens out as the incident neutron 
energy increases. Above 6 MeV competing carbon cross sections lead to an 
enhancement of pulse heights below ~ 0.5 MeVee. 



5.2 Integral detection efficiency 



The neutron yield Nn(Tn) was obtained for detection thresholds of 0.3, 0.5, 
0.7, and 1.0 MeVee from the TOF distribution (Fig. 7), converted to T„ assum- 
ing a path length of 1.5 m. Only events with T„ >1 MeV (TOF below 110 ns) 
where a software pulse-height threshold could be set above the uncertain hard- 
ware threshold range were analysed. Substituting A^„(T„) for (T„) in Eq.3 
delivered the absolute detection efficiency, but the efficiency could have been 
calculated alternatively by integration of the pulse-height spectra. 

The average efficiencies, excluding the bottom detectors of the array, are pre- 
sented in Fig. 9 and compared to the MC predictions. At neutron energies very 
close to threshold the simulation under predicts the data, especially for the 
distributions at 0.3 and 0.5 MeVee threshold. The sharp rise in the efficiency 
just above threshold is well described by the MC simulation, although here 
the data fall slightly below the model. Generally, the agreement is best in the 
region between 5-7 MeV where the efficiency function flattens out and here 
the relative difference, defined as {Sexp — £mc) I^mc^ is within 5%. At higher 
energies the statistical precision of the data deteriorates due to the rapid fall 
in fission yield. 
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Figure 8. Differential detection efficiencies (Eq.3) at neutron energies 2.25 - 
7.25 MeV, for 0.5 MeV-wide bins of neutron energy. The full-line histogram is de- 
rived from the MC simulation, while the data points are derived from the measured 
pulse-height response. Error bars denote statistical uncertainties. 

5.3 The deuteron photodisintegration measurement 

As a test of the ability of the MC model to correct for systematic scattering 
and absorption effects in neutron photoproduction, the (7, n) p cross sec- 
tion was measured with tagged photons using a D2O target. The 9 cm long, 
cylindrical target had a diameter of 6.6 cm and a significant fraction of the pro- 
duced neutrons interacted in the target volume. The Nordball detectors were 
positioned as in the ^^^Cf efficiency measurement (Fig. 1) with the target cen- 
tered at the position previously occupied by the fission source. A 18 mm thick 
plastic scintillator placed between target and detectors aided identification of 
charged particles, mainly generated by Compton scattering and pair produc- 
tion in the target. Monoenergetic photons in the range = 14 — 18 MeV were 
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Figure 9. Measured and simulated absolute detection efficiencies for the Nordball 
array. Thresholds are 0.3, 0.5, 0.7 and 1.0 MeVee. The full-line histogram represents 
the MC simulation. Error bars on the data points denote the statistical uncertainty. 

produced by the tagged bremsstrahlung technique using the MAX-lab tagging 
spectrometer [2]. The width of a tagged-photon energy bin was 340 keV and 
the rate in each plastic-scintillator element of the focal plane was ~ 400 kHz. 

The 6-8 MeV neutrons produced by two-body photodisintegration interacted 
in the detectors and bulk materials in the target-detector locality. These in- 
teractions were evaluated with the coupled GEANT3-STANT0N simulation, 
described in Sec. 4.2, to provide correction factors for the measured neutron 
yields. 

Fig. 10 shows the comparison between the present data, a high precision mea- 
surement at proton Oqm = 90° [26], and a fit to available low energy (7, p) n 
data [27]. This fit gives a very good account of the low energy data set and 
agrees well with a subsequent low energy evaluation [28]. The standard rep- 
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resentation of two-body deuteron photodisintegration is with respect to the 
proton CM angle. The good agreement between data and evaluation supports 
the ability of the simulation to take scattering and absorption effects into ac- 
count, bearing in mind that ~ 30% of the neutrons interacted in the target 
volume and also a significant fraction in the plastic-scintillator sheet. Note 
that the simple analytical approach to estimate neutron absorption in each 
material, on the basis of the reaction cross section and an exponential at- 
tenuation function, underestimates neutron losses by more than 30%. This 
approach fails to account for elastic scattering effects which are important at 
low neutron energies. 

The average cross section measured at the 5 angles was 81.4 ± 4.5(stat) ± 
9.7{sys)fih/sT and the parameterisation of [27] gives a value of 84.7±1.7/ib/sr 
where the uncertainty is due to the uncertainties in the coefficients of the Leg- 
endre expansion of the angular distribution. For the present data, the sources 
of systematic uncertainty were detector solid angle (0.7%), neutron scatter- 
ing and absorption (3.2%), detection efficiency (5%), and random-background 
subtraction (7.6%). The uncertainty for scattering and absorption was esti- 
mated by comparing neutron attenuation coefficients from other experiments 
(Sec. 5.3.1) with the predictions of the present MC simulation. 

The present measurement at neutron 6'iab = 90° is also consistent with the 
datum of [26] (Fig. 10). 



5.3.1 Application to other detector systems 

The present simulation gives an excellent account of the tagged (by the pion 
capture reaction 7r~p^n7) neutron measurement of Ref. [29]. This reported 
not only the absolute efficiency of a 16 (4x4) cell array of liquid scintillators (of 
somewhat smaller volume than the Nordball elements), but also the magnitude 
of "cross-talk" effects where neutrons scatter from one cell to another. The 
pion capture reaction produced mono-energetic neutrons of 8.9 MeV kinetic 
energy which were tagged by detection of the 129 MeV photon in a large NaT 
detector. With two-body kinematics the angular range of the tagged neutron 
beam was determined by that of the NaT detector, chosen such that the neutron 
beam directly illuminated the four central neutron cells only. The measured 
detection efficiency was 21.3±0.1(stat)±0.6(sys)% for the full array while the 
simulated efficiency with the present code was 21.6±0.1(stat)%. Magnitudes 
and shapes of the pulse-height distributions for cross-talk events were also 
well reproduced. Further comparisons of simulation predictions with neutron 
measurements [30,31] also show good agreement. 
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Figure 10. Comparison of the present ^H(7,n)p differential cross section (solid cir- 
cles), converted to the equivalent proton CM angle, for = 14 — 18 MeV with the 
evaluation of [27] (solid line) and the datum of [26] (solid square). The error bars 
on the present data denote the statistical uncertainties while those of [26] show the 
total uncertainty. 

6 Conclusions 

The scintillation response and detection efficiency of the Nordball liquid scin- 
tillator array for neutrons was measured using the well known fission, fast- 
neutron yield from a ^^^Cf source. Using the measured recoil proton pulse- 
height response as input to a MC computer model of Nordball, good agree- 
ment between simulated and measured detection efficiencies was achieved for 
neutron energies up to 8 MeV. Above this the fission yield is too low for 
statistically significant calibrations. 

The MC simulation is intended to provide systematic correction factors for 
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neutron photoproduction cross sections and this was checked by measuring 
the well known ^H(7, n) p cross section in the energy range = 14 — 18 MeV. 
Comparison with parameterisations of the 7 + d ^ p + n differential cross 
section and previous data suggest that the uncertainty of the MC calculation 
is within 5%. 

In the future we plan to convert of the present GEANT-3 based code to the 
more modern object-oriented GEANT-4 framework. This is now used exten- 
sively to model ionising-radiation interactions in matter at all energies and 
has branched into many applications including medical physics. 
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